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Abstract 

We study the muon anomalous magnetic moment = (5^ — 2)/2 in a supersymmetric U{iy 
model. The neutralino sector has extra components from the superpartners of the f7(l)' gauge 
boson and the extra Higgs singlets that break the symmetry. The theoretical maximum 

bound on the lightest neutralino mass is much smaller than that of the Minimal Supersymmetric 
Standard Model (MSSM) because of the mixing pattern of the extra components. In a model 
where the symmetry is broken by a secluded sector (the S-model), tan f5 is required to be ^ 3 

to have realistic electroweak symmetry breaking. These facts suggest that the prediction may 
be meaningfully different from that of the MSSM. We evaluate and compare the muon anomalous 
magnetic moment in this model and the MSSM and discuss the constraints on tan/3 and relevant 
soft breaking terms. There are regions of the parameter space that can explain the experimental 
deviation of from the Standard Model calculation and yield an acceptable cold dark matter relic 
density without conflict with collider experimental constraints. 
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I. INTRODUCTION 



The anomalous magnetic moment of the muon = {g — 2)^/2 is one of the most precisely 
measured physical quantities. Its current value from the Brookhaven National Laboratory 
E821 experiment is |^ 

a^(exp) = (11 659 208 ± 6) x 10^^°, (1) 

which is a 2Aa deviation from the Standard Model (SM) prediction 

Aa^ = a/,(exp) - a^(SM) = (23.9 ± 10.0) x 10"^° (2) 

when the hadronic vacuum polarization information is taken directly from the annihilation 
of e^e~ to hadrons 3] measured at CMD-2 |^]. The uncertainties involved in Eq. Q are 
7.2 X lO"^*' from the leading-order hadronic contribution [5j, 3.5 x 10~^° from the hadronic 
light-by-light scattering Q, and 6 x 10^^'^ from the experiment. The indirect hadronic 
information from the hadronic r decay gives a higher SM value that does not indicate a 
significant discrepancy with the SM (only a 0.9ct deviation)^. Recently released KLOE 
data [8] show an overall agreement with the CMD-2 data confirming that there is a 
discrepancy between the hadronic contributions from the e^e~ data and the r data obtained 
from ALEPH, CLEO and OPAL ||. 

New physics is expected to exist at the TeV-scale to resolve various theoretical problems, 
including Higgs mass stabilization, and new physics could give a significant contribution to 

to explain the above deviation lOj. There have been extensive studies of in super- 
symmetric (SUSY) models llj, which show that supersymmetry can naturally explain the 
deviation of Eq. 

The data constrains the SUSY parameters, including the sign of /i 0] and upper 
limits on relevant scalar and fermion superpartner masses [l3| . In the minimal supergravity 
model (mSUGRA) or the Minimal Supersymmetric Standard Model (MSSM), the dominant 
additional contribution to comes from the first-order radiative corrections of the chargino- 
sneutrino and the neutralino-smuon loops; it is 

A. rqnqv^ io.. m-io tan sign (/i) 
Aa,(SUSY) ~ 13 X 10 (^^^^^/^qq GeV)^ 

in the limit that all the supersymmetric masses are degenerate at MgusY 3- The 2-loop 



corrections involve sfermion subloops or chargino/neutralino subloops and are at about the 
few percent level, although full calculations are not yet complete [l5[- The discrepancy in 
Eq. (0) shows that a supersymmetry solution can be found if sign(/i) > and MgusY ^ 700 
GeV for tan/5 ^ 50, in the limit that supersymmetric masses are degenerate. The deviation 
of similarly gives constraints on the parameters of other new physics models including 
the mass of a second generation leptoquark the mass of the heavy photon in the little 
Higgs model and the compactification scale of an extra dimension [l8| . 

Given that has been a powerful tool for constraining the new physics models, due to the 
accuracy of its measurement and the SM evaluation, it is interesting to pursue what can 
tell about recently emerging models. The recent idea of split supersymmetry assumes large 



^ For a recent review of the various SM predictions and the a^t discrepancies, see Ref. 



2 



masses (e.g., 10^° GeV) for scalar superpartners (sleptons, squarks) while keeping fermionic 
superpartners (gauginos, higgsinos) at the TeV-scale |19|. The large masses of the smuon 
and sneutrino would make the chargino-sneutrino and neutralino-smuon loop contributions 
to negligible; the split supersymmetry model would be rejected if the deviation of is 
in fact real. 

Another interesting TeV-scale new physics model is the supersymmetric model [ic 



2l|. It has a structure similar to the MSSM but has an extra U{1) gauge symmetry {U {!)'), 
which is spontaneously broken at the TeV-scale by one or multiple Higgs singlets. This 
model can provide natural solutions to some of the difficulties the MSSM faces, including 
the explanation of the electroweak scale of the /i parameter (//-problem ji^]) and the lack^ of 
a sufficiently strong first-order phase transition for electroweak baryogenesis (EWBG) jj^ . 
The Next-to-Minimal Supersymmetric Standard Model (NMSSM) [2_5] can also resolve the 
//-problem but its discrete Z3 symmetry invokes a cosmolo gica l domain wall problem ji^; a 
variant which avoids this problem is discussed in Ref. 1271128^. 

Besides the bottom-up reasons to introduce an additional U{1) symmetry to supplement 
the MSSM, many new physics models, including grand unified theories (GUTs), extra di- 
mensions ji^, superstrings js^l, little Higgs [31], dynamical symmetry breaking [s^ and 
Stueckelberg mechanism models js^ predict extra U{1) symmetries or gauge bosons. The 
newly introduced particles such as the U{1)' gauge boson {Z') and the U{1)' breaking Higgs 
singlet (S) and their superpartners Z'-ino (Z') and singlino (S), alter the Higgs and neu- 
tralino spectra. The modified Higgs spectrum js^l and the neutralino relic density js^ have 
been recently studied in the U{iy model with a secluded f/(l)' symmetry breaking sector 



21|, and the difference in the predictions from the MSSM detailed. 

There have been studies of the muon anomalous magnetic moment in models with ad- 
ditional gauge groups or Eq GUT [s^, 3^. They mostly concentrated on the loops 
including the Z' or the exotic quarks in the model and their superpartners, and constraints 
were obtained on their masses or couplings. To explain the deviation, the Z' masses 
should be typically smaller than the experimental limits of Mz' ^ 500 ~ 800 GeV from 
direct searches at the Tevatron |39l. li^. 

In this paper we quantitatively study a supersymmetric f/(l)' model with a secluded 
f/(l)' symmetry breaking sector (the S-model) |2l| to see how the extended neutralino 
sector contribution to is different from the MSSM prediction. The superpotential in this 
model is 

W = hsSHiH2 + A.SiSaSs, (4) 

where hs and are dimensionless parameters. The /i-problem is solved by replacing the /t 
term by an effective /i parameter 

/ioff = hss/ \Pi (5) 

where s/-\/2 is the vacuum expectation value (VEV) of the Higgs singlet S that acquires 
a VEV at the electroweak or TeV scale. The Z' has a large mass generated by the Higgs 
singlet fields 5*1,2,3, which acquire large (TeV scale) VEVs for small because of an almost 
F and D fiat direction. The extra Higgs singlets allow /iefr to be at the electroweak scale 
while keeping the Z' heavier than the experimental limit. 



^ The required strong first-order phase transition for EWBG is allowed in the MSSM only if the light Higgs 
mass is only slightly above the LEP experimental bound and the light stop mass is smaller than the top 



mass 
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The electroweak symmetry breaking is driven by electroweak scale trilinear soft terms, 
leading to small values for tan/3 = V2/V1 (tan/3 ~ 1 to 3), while solutions without unwanted 
global minima at {H^) = typically have (S) <^ 1.5 (Hf) A small tan/3 implies a 

problem in the MSSM; for example, the light Higgs mass needs large tan f3 (for reasonable 
superpartner masses) to satisfy the LEP bound of > 115 GeV. The mixing of the Higgs 
doublets and singlets in the S'-model, however, lowers the LEP rrih bound significantly j3J|. 
Furthermore, the maximum theoretical value for rru for a given tan f3 is increased, both 
by F-terms (similar to effects in the NMSSM ji^ [s^) and by D-terms 41 1. As a result, 
tan /3 ~ 1 to 3 is experimentally allowed in the 5'-model. 

The neutralinos in the S'-model have extra components of the Z'-ino and singlinos. The 
lightest neutralino of the t/(l)' model is often very light because of the mixing pattern of 
the extra components. Eq. Q suggests that in this model may be significantly different 
from that of the MSSM. We numerically investigate the differences and obtain constraints 
on tan (3 and the relevant soft breaking terms. We do not include Z' loops since the Z' mass 
bounds from direct searches at colliders imply negligible effects. 

In Section IH] we describe the formalism to compute the chargino and neutralino con- 
tributions in the S'-model. In Section IIIH we give the numerical analysis of and make 
comparisons to the MSSM results, before the conclusion in Section HVl 



II. SUPERSYMMETRIC CONTRIBUTIONS TO WITH U{1)' SYMMETRY 



In this section we describe the neutralino and chargino contributions to in the S'-model. 
The formalism is similar to the MSSM with a straightforward extension. We assume that 
the VEVs of the S'1,2,3 are large compared to the VEVs of other Higgs fields (if°, H2 and 
S) and that their singlino components essentially decouple 21. 33|. 



Neutralino Contribution 



Ignoring the S'1^2,3 singlinos, the neutralino mass matrix in the basis W^^ H^, H2, S, 
Z'l is given by 



V 



Ml 


91V2/2 







M2 

g2Vi/2 






-giVi/2 
92V1/2 


-hss/V2 
-hsV2/\/2 



-g2V2/2 

-hsslV2 


-hsVi/\/2 



9z'Q'{HP 







-hsV2/V2 
-hsVi/V2 


9z'Q'{S)s 



\ 



gz'Q'{m)v^ 
gz'Q\m>2 

9z'Q\S)s 
Mv ) 



(6) 



where e = g\ cos dw = g2 sin 9w', gz' is the f/(l)' gauge couphng constant, for which we take 
the GUT motivated value of gz' = ^j5/3gi. Q' is the 11(1)' charge, and the anomaly-free 
charge assignments based on Eq GUT can be found in Ref. |4^ . 

246 GeV. The diago- 



The VEVs of the Higgs doublets are (if ^ ) 



V2 



with 



nalization of the mass matrix can be accomplished using a unitary matrix A^, 



A^^M^o A^ = Diag (M^o , M^o , M^o , M^o , M^o , M^o ) . 



(7) 
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FIG. 1: Supersymmetric contributions to involving charginos and neutralinos 

The first 4 x 4 (5 x 5) submatrix of Eq. © is tfie MSSM (NMSSM) limit. Due to tlie singlino 
addition, tfiere exists a kind of see-saw meclianism that makes the lightest neutralino very 
light HEi, less than 100 GeV in the case of Mi/ 3> Mi (where the Z' practically decouples 



and the mass matrix becomes the NMSSM limit) [35[. We will consider both this limit and 

that in which the gaugino mass unification relation, My = Mi = |^M2 ~ O.5M2, is 

f 2 

satisfied. 

The smuon mass-squared matrix is given by 



MIl = + + (^3m - sin^ ^h/)M| cos 2/? (9) 

Mrr = m\ + ml + (Q^ sin^ ^v^)M| cos 2/? (10) 

Mlji = m^(A;-/iefftan/5) (11) 

ML = {Ml^y (12) 
Its diagonalization can be accomplished through the unitary matrix D as 

D^MlD = Diag{Ml,Ml). (13) 

The LEP2 SUSY Working Group analysis found m^^ 95 GeV from jl nxi 2 searches 
01 . Since is small compared to supersymmetric parameters, the off-diagonal terms of 
Eq. (jlll) and Eq. ()12|) are small and hence the mixing is small. For tan/? ^ 1, both m|, 
^% ^ (95 GeV)^ are required to give m^^ ^ 95 GeV from the tree- level mass matrix of 
Eq. (jSJ, while for tan/? ^ 1, somewhat larger values of these parameters are required. We 
require both m| and m|, to be larger than (100 GeV)^ and take = in our numerical 
analysis. 

The neutralino contribution to is then |1 



where 



= aUx") + aUx") (14) 



^(X") -EE ^Re[L,,W,,]F,(§t) (15) 

j=l k=l °^ ^^^X° 
6 2 ™2 /\^2 



]=lk=l X Xi 
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with the following chiral coupling: 

= ^ {giY,,N*^ - g^N;^ + gz'Q'{^^L)N;^ D,^ + ^^Nl^^k (17) 

1 \/2m 
= -= {giY^^N^, + gz'Q'{l^R)N^,) D^u + ^N^jD^^- (18) 

= —1, Y^j^ = 2 are hypercharges and the terms with gz' coupling are the additional 
contributions from the f/(l)'. Since is small compared to the supersymmetric masses, 
we can approximate it as zero in the loop integral functions to obtain 

^i(^) = ^ ^,J l-x^ + 2x\nx) (19) 

F2{x) = ^^^J -x^ + 6x^-3x-2-6xlnx). (20) 



B. Chargino Contribution 

The chargino mass matrix is given by 

M,. = ( , ^f'"-'""). (21) 

^ \\/2Mw cosp hsS/^/2 J ^ ' 

It is essentially the same as in the MSSM except that /i is replaced by /igff = hss/\/2. M^± 
can be diagonalized by two unitary matrices U and V as 

U*M^±V-^ = Dtag{M^±,M^±). (22) 

The LEP light chargino mass limit M^- 104 GeV 0] gives constraints on M2 and fics 
for a fixed value of tan f3. The sneutrino mass-squared is given by 

= m| + Tg^Ml cos 2/5. (23) 

As in the MSSM calculations, we do not include the right-handed neutrino or its superpart- 
ner. 

The chargino loop contribution to is 

^mIx') =ai(x") + aj(x") (24) 

where 

alix-) = E E ^^^^[^^-^.*^]^3(^) (25) 

j=l k=l Xj Xj 

2 1 ™2 Tur2 

4(x') = - E E + Mj^) (26) 

j=l k=l Xj 



with the chiral ^-i^n-X couplings 



V2', 



Ln = "^U;,, R,, = -g,V,,. (27) 
f 1 
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FIG. 2: (a) Aa^, in the MSSM [filled circles] and the NMSSM hmit (Mi/ » Mg) of the 5-model 
[dark shading] for tan /? = 2.5 and rngmuon = 100 GeV. For the 5-model, hg = 0.75 and ?7-model 
charges are assumed. The region outside the NMSSM region has Aa^ x 10^° < 13.9, while the 
island inside the NMSSM region (around M2 ~ 100 GeV) has Aa^ x 10^° > 33.9. The hghtly 
shaded area outlined by solid curves is excluded by the LEP chargino mass limit M^- > 104 GeV. 
The Attf, have values 13.9 < Aa^ x 10^° < 33.9 favored by Eq. @. The models give similar 
results and all the models allow small tan 3 for rrismuon small (~ 100 GeV). (b) The solution space 
for the acceptable neutralino relic density |35| and the {g — 2)^ deviation are shown together for the 
same NMSSM hmit. The filled squares have the WMAP 3a allowed range of 0.09 < O^o/i^ < 0.15, 
the open circles have fi^o/i^ < 0.09 and the crosses have 0.15 < Q^oh"^ < 1.0. 

The loop integral functions are 

F^ix) = -l-^-—{3x^-4x + l-2xHnx) (28) 
/ (a; — l)'^ 

FAx) = --- — ^—-{2x^ + 3x^-6x+l-6x'^\nx). (29) 
6 [x — 1)^ 

III. ANALYSIS 

The MSSM can explain the 2Aa deviation between the E821 experiment and the SM 
prediction for most values of tan f3. In this section, we compare predictions from the MSSM 
with those of the S-model. 

Figure|21(a) shows Aa^ in the M2-S (also M2-/i) plane for the MSSM [filled circles], the S- 
model and its NMSSM limit [dark shading], with tan/? = 2.5 and rngmuon = ""^l = ''^^ = 100 
GeV. The comparison is made for /i = /ieff. The corresponding plot for the S'-model with 
Mil = Ml is hardly distinguishable from that of the NMSSM limit so it is not shown in the 
figure. The plots are the parameter spaces that can give values Aa^ x 10^° = 13.9 ~ 33.9 
favored by Eq. Q. The parameter space excluded by the LEP chargino mass limit of 
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FIG. 3: Maximum Aa^ versus (a) tan/3 and (b) Tn-smuon in the S'-model with Mi/ = Mi. The 
MSSM result and the NMSSM hmit {Mr > M2) are nearly indistinguishable from the Mi> = Mi 
curve and are not plotted. M2 is scanned from 100 GeV to 1000 GeV, and fi{fics) from 100 GeV 
to 1000 GeV. 

'T^smuon is Scanned from 100 GeV to 500 GeV for (a), and tan/3 from 1 to 3 for (b). 
Horizontal dashed-dot lines are boundaries of the measured deviation Aa^ x lO^'^ = 13.9 ~ 33.9 of 
Eq. ©. 

M^- > 104 GeV is shown as the lightly shaded area outlined by solid curves. Throughout 
the following analysis (Figures IHl and we do not include the parameter points that violate 
the LEP chargino mass constraint. This figure shows that both the MSSM^ and the S'-model 
(including the NMSSM limit) can explain the Aa^ data with small tan (3 while satisfying 
^smuon ^ 100 GeV from the experimental bound on the scalar muon mass. The S'-model 
gives only a slightly larger area in the parameter space than the MSSM. For the S'-model, 
hs = 0.75, and the Eq motivated r^-model'^ charge assignments are assumed |4^. 

It is interesting to note that in the Mir 3> M2 case^ (NMSSM limit), a sizable part of 
the 2.4(7 deviation solution area overlaps the solution area which reproduces the observed 
cold dark matter (CDM) relic density in the same framework and limit js^ as shown in 
Figure H (b). The filled squares have the WMAP 3cr allowed range of 0.09 < Q^oh'^ < 0.15, 
the open circles have Vt^oh? < 0.09 and the crosses have 0.15 < Vt^oh'^ < 1.0. In the event 
that the Aa^ deviation from the SM is real, this CDM result enhances the viability of the 
model. The chargino mass allowed by both the deviation and the relic density is 104 
GeV ^ <^ 220 GeV, where the lower bound is the present experimental limit. Since 
the chargino mass could be only slightly larger than the present LEP limit, a search for the 



^ It should be emphasized that the MSSM cannot accept this small tan f3 if we consider other constraints: 
the light Higgs mass would be too small to be compatible with the LEP Higgs mass bound of m/i > 115 
GeV. 

The ry-model is a U{iy model that is produced with a unique set of charge assignments when Eq is broken 
to a rank- 5 group. 
^ The CDM expectations for smaller My have not yet been examined. 
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SUSY trilepton signal at the Tevatron Run II will be very interesting j46[. It should be 
noted that only the Z-pole annihilation channel was considered to show this model could 
reproduce the acceptable relic density. There could be a larger solution space when more 
channels are considered. 

We now consider the limits on tan/? and rngmuon that allow the favored range of Aa^. 
Figureini(a) shows the maximum Aa^ as a function of tan (3 in the S'-model (with My = Mi). 
The MSSM curve is nearly the same, as is the NMSSM limit (Mi/ ^ M2) of the S'-model. 
M2 is scanned from 100 GeV to 1000 GeV, /i(/ieff) from 100 GeV to 1000 GeV (negative fx is 
nearly irrelevant for producing positive Aa^), and msmuon from 100 GeV to 500 GeV. This 
figure demonstrates the almost linear dependence on tan/5, as in Eq. 0, and shows that 
even with small tan/3, both models are able to produce the favored values of Aa^. Smaller 
tan j3, however, has less supersymmetric parameter space that explains the deviation. The 
maximum values of Aa^ have msmuon ~ 100 GeV (the lowest scan value). 

Since the chargino contribution is the same in the MSSM and the S'-model, the dominant 
difference in would come from the lightest neutralino (x?) contribution. The mass and 
the coupling of Xi the S-model is similar to that of the MSSM. In the NMSSM limit 
(Ml/ 3> Ml), Xi is often significantly lighter than that of the MSSM but it is mostly 
singlino-like which couples to the muon only through mixing; X2 is mostly similar to Xi of 
the MSSM. 

Figureini(b) shows the maximum Aa^ as a function of msmuon in the S'-model; the MSSM 
curve is nearly identical when tan f3 is scanned only from 1 to 3. The maximum points have 
tan/3 ~ 3 (the highest scan value). A large value of tan/3, as more general U{iy models 
allow, would always ensure a sufficiently large maximum Aa^ for a given msmuon- This 
figure demonstrates the roughly inverse- squared dependence on rrismuon as in Eq. Q. It also 
shows that, for small tan/3, the scalar muon should be light (100 GeV ^ M^ ^ 180 GeV) 
to explain the deviation. This could in principle lead to a concern for the neutralino 
cold dark matter candidate, since the scalar muon could be light enough to be the lightest 
supersymmetric particle (LSP). A charged LSP would confiict with the observational absence 
of exotic isotopes. In the S'-model, however, the lightest neutralino is usually very light, e.g., 
M^o <^ 100 GeV, while it produces the acceptable range of the relic density for a wide range 
of the parameter space j3^. One can therefore have a light slepton with an even lighter 
neutralino as the LSP. A light slepton will be observable in future accelerators such as the 
CERN Large Hadron Collider (LHC) and International Linear Collider (ILC). Light sleptons 
could be detected easily at a linear collider of a moderate energy of 500 GeV. 

Figure 0] shows the maximum Aa^ as a function of M2 in both models. For the dotted 
curves, we do not impose the gaugino mass unification assumption of Mi/ = Mi ^ O.5M2, 
and we scan Mi from 50 GeV to 1000 GeV. This figure shows that independent Mi can 
make quite a difference for a fixed M2 when M2 is large. For Figure El a relaxation of 
gaugino mass unification would not make much difference since the maximum Aa^ mostly 
happens for small M2^. In the S'-model we do not relax Mi/ as a free parameter for the 
dotted curve, but rather take the NMSSM limit (Mi/ ^ M2) and relax only Mi. Unless 
Mil is very small (smaller than the Mi scan limit of 50 GeV), it would not increase the 
maximum Aa^ significantly. For a wide range of M2 (and practically for any M2 in the case 
that Ml is a free parameter), both models can produce the favored Aa^. The maximum 
Aa^ have tan/3 ~ 3, msmuon ~ 100 GeV (and Mi ~ 50 GeV for the dotted curves). 



^ Small AI2 results in a small chargino mass in Eq. (|21|l . and large Aa^ in Eq. H24|l . 
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FIG. 4: Maximum Aa^ versus M2 in (a) the MSSM and (b) the S'-model. The dotted curves 
take Ml as a free parameter while the sohd ones fohow the gaugino mass unification relation 
[My = Ml ~ O.5M2). ^(^eff) is scanned from 100 GeV to 1000 GeV, Mi (for dotted curves) from 
50 GeV to 1000 GeV, tan /3 from 1 to 3, and m^mnon from 100 GeV to 500 GeV. When Mi is a free 
parameter, we take My ^ M2 (NMSSM limit) for the S'-model (dotted curve). In both models, 
the maximum values of Aa^ occur for Mi ~ 50 GeV when Mi is taken as a free parameter. 

IV. CONCLUSION 

Unlike the MSSM, the electroweak symmetry breaking condition in the S'-model requires 
small tan/3. Moreover, the LEP smuon mass bound requires slepton masses above about 
100 GeV. On the other hand, the 2.4(T deviation of the muon anomalous magnetic moment 
favors large tan (3 and small slepton soft terms. The Aa^ determination gives the most severe 
constraint on tan/5 and slepton masses in the S-model; nonetheless, the S'-model can still 
explain the deviation while satisfying the chargino and smuon mass limits. Since the mass 
of the light smuon is constrained to be less than 180 GeV, it would be easily observable at 
the next generation colliders. It is remarkable that the parameter space that explains the 
deviation has a sizable overlap with the parameter space that produces an acceptable cold 
dark matter relic density even when only the Z-pole annihilation channel is considered. The 
common solution space implies a lighter chargino upper mass bound of 220 GeV, though the 
common solution space would increase when annihilation channels in addition to the Z-pole 
are included for the neutralino relic density calculation. More general t/(l)' models that do 
not require small tan/? could explain the deviation in a wider range of the parameter 
space. 

Even though the lightest neutralino in the S'-model is often very light, the difference of 
the Aa^ predictions of the MSSM and the f/(l)' models for comparable parameters is not 
large. This is because the lightest neutralino is mostly singlino-like when it is lighter than 
that of the MSSM, and it couples to the muon only through mixing; the properties of the 
other neutralinos are quite similar to those of the MSSM. The relaxation of gaugino mass 
unification makes a sizable difference in both models. 
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The contribution of the Z'-\oop is suppressed by the large Z' mass, which is constrained 
by the CDF hmit of Mz' ^ 500 ~ 800 GeV, with the hmit depending on the model 0|. In 
the case that the right-handed neutrinos are f/(l)' charged and form Dirac particles, more 
severe constraints of Mz' ^ multi-TeV are deduced from Big Bang Nucleosynthesis (BBN) 
|47l |. Other possibilities for neutrino mass in these models are discussed in Ref. jisj ]. 
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